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Abstract 

The local zero structure of a smooth map may qualitatively change, when the map is sub¬ 
jected to small perturbations. The changes may include births and/or deaths of zeros. The 
qualitative properties are defined as the invariances of an appropriate equivalence relation. The 
occurrence of a qualitative change in the zero structures is called a bifurcation and the map is 
named a singularity. The local bifurcation analysis of singularities has been extensively studied 
in singularity theory and many powerful algebraic tools have been developed for their study. 
However, there does not exist any available symbolic computer-library for this purpose. We 
suitably generalize some powerful tools from algebraic geometry for correct implementation 
of the results from singularity theory. We provide some required criteria along with rigorous 
proofs for efficient and cognitive computer-implementation. We have accordingly developed 
a Maple end-user friendly library, named “Singularity”, for an efficient and complete lo¬ 
cal bifurcation analysis of real zeros of scalar smooth maps. We have further constructed a 
built-in help for Singularity and provided a comprehensive user-guide. Recent progresses 
in Maple have been used so that Singularity generates the list of inequivalent persistent 
bifurcation diagrams for parametric singularities. The approach is independent of the singular¬ 
ity's codimension (a measure of its complexity). The main features of Singularity are briefly 
illustrated along with a few examples. Our Maple library is not only useful for research goals 
and engineering applications involving singularities but also for pedagogical purposes. 

Keywords : Singularity and bifurcation theory; Ideal membership problem; Standard and Grobner 
bases. 

2010 Mathematics Subject Classification: 37G10; 13P10; 58K50; 58K60. 

Public access: Singularity will be released for public access once this draft paper is peer-reviewed 
in a refereed journal. 


Corresponding author. Phone: (98-31) 33913634; Fax: (98-31) 33912602; Email: mgazor@cc.iut.ac.ir. 




M. Gazor and M. Kazemi 


Bifurcation analysis for local zeros of smooth maps 


2 


Many real life problems may result in the analysis of local zeros (around a zero solution, named 
a base point) of a smooth map 

/:R n xr^ R", f(x, a) = 0, for x G M n , a G M m (0.1) 

We refer to x by state variables, n by state dimension and a by parameters. Note that locating 
singular base points of a smooth map is related to finding roots of nonlinear systems and is not the 
purpose of our work here; see HZIEHEII. Hence, we may assume that the base point is the origin 
and /(0, 0) = 0. Equation (10.11) may demonstrate a surprising change on the solution set when the 
parameters vary. This occurs when the Jacobian matrix of / does not have a full rank. In this case 
we say that / is singular. 

Equation (10. 1 1) may appear by direct mathematical modeling of a singular engineering problem 
or indirect, through reduction methods such as Liapunov-Schmidt reduction ; e.g., see [30], [31, Pages 
156-162] and [29] Chapter VII]. For example, Equation (10.11) appears in the study of equilibria 
and limit cycles of dynamical systems or steady-state solutions of partial differential equations. 
In fact, the theory described here is known as a “natural framework” for equilibrium bifurcation 
theory ; see [30]. Using Liapunov-Schmidt reduction, we can reduce the state dimension so that the 
Jacobian matrix at the origin is the zero matrix. Thus in this paper, we assume that 

df 

n — 1 and —— (0,0) = 0. (0.2) 

ox 

We will deal with the case of multi-state dimensional problems in a future project. When / is a 
singular map and the parameters a vary, the number of solutions for Equation (10.11) may change 
and any of such changes is called a bifurcation. The equation f(x,a ) = 0 is called a bifurcation 
problem and the set 

{(x,a) : f(x,a) = 0} (0.3) 

is called a bifurcation diagram. 

Singularities (including those different from here) and the theories to deal with them have 
appeared in several mathematical disciplines with the same or similar keywords. These theories have 
a long history dating back to the original works of Hassler Whitney, and then John Mather, Rene 
Thom and V.I. Arnold in late 1960s. Depending on the context, space and equivalence relation, they 
are faced with essentially different tools, algebraic and geometric structures, applications and/or 
objectives; e.g., see [2ll7][l5 ] l29 ][32]HTl H5] . For an instance a geometric singularity refers to a variety 
whose tangent space is not regularly defined; e.g., a self intersecting curve or a cusp; see [I51I32]. 
This leads to a theory of its own and it is named singularity theory. Arnold in HH uses pt- and right- 
equivalences for his normal form classification of smooth functions. Strategy equivalence is used for 
the analysis of strategy functions in [T5]. Right-left equivalence is applied in [7] to treat bifurcations 
of Hamiltonian systems. The other instance is catastrophe theory [44] and usually uses almost 
identical keywords as ours. Catastrophe theory addresses basically different qualitative properties 
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using right equivalence than the qualitative properties we pursue here using contact equivalence. All 
of these mentioned subject areas follow different objectives and goals than what ours do. In this 
paper we only refer to the methodologies in dealing with the local bifurcation analysis of real zeros 
of smooth maps (10.111 as singularity theory; see pj^l [l8H20l[29llMll38l[391ITT] for our main references 
of this subject. 

Many powerful algebraic tools have been developed for local bifurcation analysis of zeros in 
Equation (10. Ill . Armbruster [3] proposed a cognitive use of Grobner basis and encouraged a system¬ 
atic implementations of the existing results in a computer. Yet to the best of our knowledge, there 
does not exist any end-user symbolic library available (neither free nor paid) to everyone for the 
local bifurcation analysis of zeros of smooth maps. This is a long overdue contribution despite its 
importance and wide applications. In the last two decades, there has been a considerable progress 
in development of computer algebra systems so that an efficient symbolic implementation of the 
results in singularity theory is now feasible; see [30,31j for numeric approaches. 

Singularity theory uses two fundamental notions of intermediate order terms and high order 
terms for normal form computation. Recall that theory of normal forms here is to simplify singular 
maps and qualitative study of their zeros. Intermediate order terms have a complicated structure 
and do not seem to hold a convenient algebraic structure. Therefore, Grobner basis type of tools 
(i.e., its local versions) alone are not sufficient to simplify removable intermediate order terms in 
a singular germ. Furthermore as far as our information is concerned, there does not yet exist a 
complete algebraic characterization for high order terms in many cases of multi-state variables. This 
does not mean that high and intermediate order terms are computationally inaccessible but it does 
mean crude uses of Grobner basis type of tools fail to address normal form computations. (Note that 
the basic and principal ideas in singularity theory start with normal form computation.) Indeed, 
the approach needs to be refined along with nonlinear normalization techniques. Our approach 
efficiently simplifies removable intermediate order terms and is adoptable for simplification of high 
order terms even for the generalization of our results to multi-state variable cases. 

Our main contributions here are some natural generalizations of known tools from algebraic 
geometry for correct symbolic implementation and providing certain criteria along with rigorous 
proofs for alternative and efficient symbolic implementation, an efficient symbolic implementation 
of the results of singularity theory (one dimensional state variable) in Maple, developing an end- 
user friendly MAPLE library named Singularity, writing a user guide ED and constructing a 
comprehensive built-in MAPLE help for Singularity. 

A contribution here is to generalize some existing techniques and concepts from algebraic ge¬ 
ometry to the context of (locally) smooth maps (germs) for a correct implementation. Due to the 
infinite nature of Taylor series of smooth maps, the computations are performed modulo a given 
degree. We provide a sufficient condition for a given degree whose truncation does not lead to 
error. The default work of Singularity tests the condition and do the computations modulo an 
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optimal degree. However, this approach adds a computational cost. Further, smooth maps involve 
flat functions (functions with zero Taylor series) and this may cause unnecessarily complicated for¬ 
mulas. Thereby, it is fundamentally helpful to use a ring smaller than the ring of smooth maps 
when it is feasible. Unlike the ring of formal power series, the associated computations in the ring 
of polynomials or fractional maps are exact and no truncation is required. We provide conditions 
along with rigorous proofs for the possible efficient implementations using ideals generated in ei¬ 
ther the rings of polynomials, fractional maps or formal power series. Singularity is adopted 
accordingly. We accordingly provide guidelines for efficient symbolic implementation of the existing 
results from singularity theory; see (3lJ Chapters 6-7], [HJ Sections 6.2 and 6.3], [29[ Chapters 1-4] 
and [3llTl [T9ll39] . This simply illustrates an interesting and realistic application of computational 
algebraic geometry to the equilibrium bifurcation analysis of real world problems. 

Maple is one of the two major (along with Mathematica) symbolic computer algebra sys¬ 
tems. Besides, Maple is user-friendly and is widely used for Mathematics and engineering research 
as well as educational purposes. Thus, Maple is a good choice for end-user libraries in general, 
and in particular for Singularity due to various applications of singularity theory in different 
engineering disciplines. Further, many useful and advanced techniques from computational alge¬ 
braic geometry such as Grobner basis, elimination ideals for polynomial ring, regular chains and 
triangular decomposition method have already been developed and are now available in built-in 
packages of Maple 18. So far, the most efficient Grobner basis computation is due to Maple 
F 4 algorithm implemented by Jean-Charles Faugere (see [33]) while regular chains and triangular 
decomposition method are only available in MAPLE. Singularity has benefited from an efficient 
use of these capacities of Maple. Besides, we have already developed some Maple programs for 
(parametric and orbital) normal form computation of singular differential equations [22l ITT] and 
their integration with Singularity shall lead to a toolbox for local bifurcation control and analysis 
of singularities. Singularity has been tested by all scalar examples (nonsymmetric and without 
modal parameters) and classifications given in [29], ISlj, 38j,[39] and a few (differences, error or not 
already reported data due to computational burden) are verified in our favor. We use identical 
notations and terminologies from [T^[29ti_31:, 32 v 38|[39| as far as it is feasible. 

Singularity computes a variety of algebraic structures associated with singular scalar maps 
including tangent and restricted tangent spaces, high order term ideals, and the intrinsic ideals 
associated with ideals of both finite and infinite codimension. Our Maple library derives low and 
intermediate order terms, normal forms, universal unfolding, and transition sets. Singularity 
efficiently simplifies intermediate and high order terms. It, further, generates persistent bifurcation 
diagrams (plot or animation), and estimates the transformations transforming contact-equivalent 
scalar maps to each other. Finally, Singularity solves the recognition problem for normal forms 
and universal unfoldings. An interesting capability of Singularity is the classification of persistent 
bifurcation diagrams by generating an automatic list. This latter capability is in fact enabled by 
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using a powerful built-in MAPLE 18 package called RegularChains [ 12] . 

In order to develop Singularity, we have discussed some modified concepts (suitable in our 
context) from computational algebraic geometry including division remainders, standard bases, 
elimination ideals, ideal membership problem and colon ideals for the local rings of fractional (germ) 
maps, formal power series and smooth maps. These are accordingly implemented in Singularity. 

The rest of this paper is organized as follows. Singularity theory and bifurcation analysis of 
Equation (10. II) - (10.211 is discussed in Section [lj We further explain how singularity theory is related 
to ideal membership problem in algebraic geometry. Section [2] describes how to treat the ideal 
membership problem. In this direction, computational algebraic tools such as standard and Grobner 
bases for ideals in three different rings, and the concept of finite codimension ideals are introduced. 
Intrinsic ideals and their associated ideal representations are discussed in Section [3] We further 
explain a procedure for computing the intrinsic part of an ideal. Section 0] gives our suggestions on 
how to implement some objects and results from singularity theory. These implementations include 
high order term ideals, tangent spaces, transition sets, persistent bifurcation diagram classifications, 
normal forms and the universal unfolding. The capabilities of the main features of Singularity 
along with a few examples are sketched in Section [5] Finally, our future and in-progress projects 
are outlined in Section [6] 

1 Introduction 

Due to the local nature of the problem ( 10.Ill -( 1(121) . we recall the notion of smooth germs around 
a base point. Two maps are considered as germ-equivalent when both maps are locally identical; 
more precisely, when there exists a neighborhood of the base point so that both maps are equal on 
the neighborhood. A germ is a germ-equivalence class of a smooth map. We denote <§ for the set of 
all scalar smooth germs whose base point is the origin. From now on we merely work with elements 
of <§ rather than a scalar smooth map; see [31, 156]. 

Following [3T], Chapter 7] and [29], we study the local zeros of maps when there is only one 
distinguished parameter denoted by A, i.e., 

j:lxK—» M, g[x, A) = 0, g(0, 0) = g x (0, 0) = 0, and m := 1. (1.1) 

The effect of additional parameters may be treated by study of their small perturbations. The 
main goal of this theory is to classify qualitative types of Equation (11.111 and its arbitrary small 
perturbations. In order to achieve this goal, we first define a qualitative property as a property that 
is invariant under an appropriate equivalence relation. Here, we use contact-equivalence relation: 

• We say that the germs g and h are contact-equivalent when there exist a smooth germ S(x, A) 
and diffeomorphic germs X and A such that 


g{x,\) = S(x, \)h(X(x, A),A(A)), 


(1.2) 
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where S(x, A) > 0, X x (x, A) > 0 and A'(A) > 0. 

Other alternative equivalences have also been reported in the literature based on their applications; 
e.g., see mm where right-left and strategy equivalences have been used. Right equivalence is also 
used to study the isolated hypersurface of geometric singularities and to study the bifurcation 
analysis in catastrophe theory, see [7|[32], 144] . The contact-equivalence is here used, since it is the 
most natural equivalence relation preserving the zero structures. 

Now we briefly describe the approach in singularity theory and how the computational algebraic 
tools may help to implement them. To study the local zeros of (II.ip . we choose a representative 
(say /) from contact-equivalent class of g that is considered to be the simplest for the analysis and 
call it a normal from. In order to compute the normal form of a singular germ, we need to compute 
certain ideals in the local ring of smooth germs. This signifies the importance of the well-known 
ideal membership problem in algebraic geometry, that is, deciding what kinds of germs belong to 
an ideal generated by a given set. One may study the zero structures of the normal forms and then, 
conclude about the solution behavior of Equation (II. ip . For instance, let 

g(x, A) := exp(x 2 ) + 2 cos(x) — 3 + sin(A). (1.3) 


Using the command Normalf orm(g, 5) in Singularity, we obtain its normal form by f(x, A) : = 
j^x 4 + A. The bifurcation diagrams of / and g are depicted in Figure 1(a) and 1(b). Here, 


Transformation(g, f , 5) provides an approximation (A + x + A 2 + Xx, 1 + |A 2 — ^A 3 — ~\ 2 x — 

|r 4 ) modulo degree 5 for (A"(r, A), S(x , A)) which 
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transforms g into / via Equation (II.21b The locally invertible transformation X(x, A) sends the 
bifurcation diagram of / into that of g. 

Real life problems can not be perfectly modeled by a system of equations and imperfections 
are always inevitable. Furthermore, the singularity of a germ g implies that the zeros of a small 
perturbation of g , say 


hfx, A, /3) = 0, where h(x, A, 0) = g(x, A), 


(1.4) 


may behave substantially different than what the zero structure associated with g = 0 does. The 
parameterized germ h in (j 1.4ft is called an unfolding of g. Hence, modeling imperfections and the 
possible existence of additional parameters in a model are the main obstacles of simply using normal 
forms of a singular germ for the qualitative understanding of a real life problem. Thus, the approach 
needs to be refined through the notion of universal unfolding. In fact, we are interested to find a 
parameterized family like 

G(x, A, a) — 0, a e M p , (1.5) 

such that for any small perturbation ep(x, A, e), the germ g + epfx, A, e) would be contact-equivalent 
to G(x, A, a(e)) for some germ a(e). We call such germ G a versal unfolding of g. A versal unfolding 
with the minimum possible number of parameters is called a universal unfolding for g\ see m 
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(a) 







Figure 1: The vertical axes stand for the state variable x. Diagrams 1(a) and l(b)| depict bifurcation 
diagrams of g in Equation () 1.3 f) and its normal form /, respectively. The diagram 1(c) shows the transition 
set and the regions 1-3 associated with Equation (11.41) . The second row illustrates the persistent bifurcation 
diagrams associated with these regions. 


Definitions 6.4.2 and 6.4.3, Page 176] and [22, Definitions 1.1 and 1.3, Pages 120-121], The number 
of parameters in a universal unfolding is named the codimension of g. The universal unfolding 
accommodates any possible modeling imperfections, any arbitrary small perturbation and also the 
existence of any possible number of parameters (in addition to the distinguished parameter A). In 
order to derive the universal unfolding of a given singularity, the computation of a vector space called 
tangent space or instead a basis for its complement is required. Using Singularity’s command 
UniversalUnf olding(g) for g in (j 1.3 [1 gives rise to 

G(x , A, a i, «2) := r 4 + A + a\x + a^x 1 . (1.6) 

Bifurcation diagram classification of the universal unfolding gives an insight to the zero structure 
of a germ and any of its perturbations. This is studied by the notion of persistence in the bifurcation 
diagrams. In fact a bifurcation diagram is called persistent when all its small perturbations remain 
self contact-equivalent, and otherwise it is called nonpersistent. Finding nonpersistent systems and 
their associated subset of the parameter space (the so-called transition set E) play a central role in 
this classification. More precisely, all parametric germs associated with parameters in a connected 
component of E c , complement of £, are contact-equivalent. Therefore by choosing a parameter from 
each connected component of E c , a complete list of persistent bifurcation diagrams modulo contact- 
equivalent is obtained. The non-persistence may either originate from a local nature or be caused 
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by the singular boundary conditions. Local nonpersistent bifurcation diagrams are determined with 
families of germs associated with three semi-algebraic parameter spaces of codimension one; i.e., 
bifurcation 0B, hysteresis and double limit point 2; see [291 Page 140] for details. Nonpersistent 
germs associated with boundary conditions add extra complications into the solution dynamics, 
when bifurcation diagrams cross the boundary; see [29, Pages 154-165]. Given our description 
for any finite codimension singularity, the connected components of the complement set of the 
transition set, i.e., E c , provide the qualitative classification of persistent bifurcation diagrams. The 
command TransitionSet(G, [x, A, aq, a 2 ]) f° r G in (1 1.6 1) gives SS := 0, := {(aq,a 2 )| 8 a 2 — 

—27al},2 := {(aq,a 2 ) | aq = 0,ai 2 < 0}, and E := 2? U 2. The transition set E is plotted 
in Figure 1(c) PersistentDiagram(G) plots a complete list of all contact-inequivalent types of 
persistent bifurcation diagrams for G; see the second row in Figure [Q 


2 Ideal membership problem and tools from algebraic ge¬ 
ometry 

Given our description in the previous section, we mainly need to address the ideal membership 
problem, in the ring of smooth germs. We call an ideal basis for the (finite) generators of a finitely 
generated ideal. The ideal membership problem refers to the question on whether or not an element 
belongs to a predefined ideal. Although a finitely generated ideal is defined by a ideal basis, yet it 
is not an easy task to understand whether a given element belongs to the ideal or not. There are 
two different ways to facilitate the ideal membership problem. One is to find a convenient ideal 
presentation from which the ideal membership can be easily understood. This is feasible for certain 
ideals and will be addressed in Section [3] The second is a computational approach based on a 
division algorithm. Dividing a given element by the ideal basis (generators), the element certainly 
belongs to the ideal when the remainder is zero. However for non-zero remainders with respect to 
an arbitrary ideal basis, the division may not help to conclude about the membership of the element 
in the ideal. B. Buchberger in 1965 introduced an ideal basis (called Grobner basis ) for ideals in the 
polynomial ring on which zero/nonzero remainders would indeed conclude the ideal membership; 
see mH- 

Armbruster and Kredel [4] suggested a cognitive use of Grobner basis for computing universal 
unfolding; also see [3]. Then, Gatermann and Lauterbach jT9] extended the tools for equivariant 
bifurcation problems. Wright and Cowel in [13] noticed that a local version of Grobner basis is 
indeed the appropriate tool to work with it in this theory. In fact, naive uses of Grobner basis 
yield wrong results for many singularities (e.g., see Remark 12.171b The main reason is that Grobner 
basis may work with polynomial germs, while the ring of smooth germs is a considerably larger ring 
than the polynomial germs. Recall that a set may generate a larger ideal in a ring than what it 
generates in a subring. Since Wright and Cowel’s remark, the only result in this direction is due 
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to Gatermann and Hosten [2D] which is a fundamental contribution. They used (mixed) standard 
basis for “mixed modules’’'’ over fractional maps and multi-dimensional state variables. They did not 
discuss computations with smooth maps and Taylor series truncations. Therefore, their suggested 
algorithms are limited to treat certain (finite codimension) fractional maps. Their algorithms are 
useful for possible implementations in a computer algebra system to simplify fractional maps with 
finite codimension. However, they are not sufficient to compute normal forms. The reason is that 
they attempted neither to simplify removable intermediate order terms nor to address the issue 
facing high order terms for multi-dimensional cases (there does not yet exist a complete algebraic 
characterizations for high order terms); for example see [20) Example 53]. 

The principal ideas of Grobner basis have been extended for localizations of the polynomial ring 
and formal power series. Here, we are concerned with local rings and the terminology “Standard 
basis” (not Grobner basis) is usually used in the context of local rings. Following Wright and Cowel’s 
remark, Gatermann and Hosten m used the localizations of the polynomial ring (a local ring) to 
circumvent the problem. The localization of the polynomial-germ ring gives rise to the fractional 
germs, i.e., germs with a fractional representative in their germ-equivalent class. Although the set of 
all fractional germs is a local ring, yet it is still a much smaller ring than the ring of smooth germs. 
The other alternative is the local germ ring of all formal power series. This is a larger ring than 
the fractional germs and it is, perhaps, suitable due to Borel lemma. (Borel lemma indicates the 
existence of a one-to-one correspondence modulo flat functions between smooth germs and formal 
power series through their Taylor series expansion.) One of the main obstacles in working with 
the local rings is termination of the algorithms. There are methods in the literature for computing 
the standard basis (for fractional maps using Mora normal form) so that they solve the problem of 
algorithm terminations. Yet for the case of formal power series, no computer program can compute 
and store infinite formal power series expansions and thus, their truncations up to certain degrees 
are mostly unavoidable. In other words, the computations are performed modulo a sufficiently high 
degree (we will justify this in Theorem 12.161 ([part 4)). In this section we investigate and compare 
the computations in the local (germ) rings of polynomials, the fractional germs and formal power 
series with those in the ring of smooth germs. We discuss the circumstances on which they can 
be alternatively used. This helps to efficiently use the different algorithms and yet ensure about 
correctness of the results. 

2.1 Standard and Grobner bases for ideals 

Let K be a held of characteristic zero; in particular we are interested in the field of real numbers. 
For our convenience, we simply identify a given germ with a convenient representative of that germ. 
For instance, we talk about the polynomial germ ring over the field K and denote it by K[x, a] while 
we mean the ring of all smooth germs whose germ-equivalent class have a polynomial representative. 
The quotient ring of <§ over the ideal of all flat germs is denoted by K[[x, a]] due to the fact that it 
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is ring-isomorphic to the ring of formal power series. Thus, we call K[[x, a]] the germ ring of formal 
power series. Further, we identify members of K[[x, a]] with the infinite Taylor expansion of their 
representative. Since K[x, A] and K[[x, A]] are Noetherian rings, we can guarantee termination of 
most algorithms during the computations. 

Denote deg(X") := ay + ... + a n when X a := x^x'f 2 ■ ■ ■ x“", X := ( 27 , x 2 , ■ ■ ■, x n ) and a := 
(ai, a 2 ,, a„.) G Z™ 0 . Any expression like X° is called a monomial germ while a term in S means 
a monomial germ along with its coefficient. We define the lexicographic ordering -<i ex on monomial 
germs X a as follows: 

X a -<i ex X' 3 when a* — j3i is negative for i := inf {j \ otj — / 3j ^ 0}. 

Definition 2.1. A local order -< is a total ordering (every two terms are comparable) and further¬ 
more, 

• For any a, /3, 7 G Z> 0 , the condition X a -< XP implies X a+1 -< X^ +1 . 

• Xi -< 1 for all i = 1 ,,n. We further assume that 0 -< X a for any a. 

Due to Dickson’s lemma (see [ 6 } Page 251] and [Hi Page 71]), every (infinite) set of monomials 
have a maximum with respect to any arbitrary local order; here, the condition X^\X a implies 
X a -< X 13 . An important example of a local order is anti-graded lexicographic ordering -< a i ex defined 
by 

( deg(A"“) > deg(X /5 ), 

^ alex X fj yy J 

[ deg(X") = deg(A^) and A Q ^ iex X&. 

The localization of the polynomial germ ring is defined as 

^:={^|/,SGA>,A],s( 0 , 0 )^o} 

whose unique maximal ideal is generated by x and A. It is common to denote & with K[x, A]( 2 , i a)- 
We call SS the ring of fractional germs. Throughout this paper, we denote 

TZ for either of the local rings K[[x, A]],^ 1 , or S, 

unless it is explicitly stated. 

Definition 2.2. Let / G TZ and -< be a local order. The infinite jet and fc-jet of / are defined by 
its Taylor series expansion around the origin and denoted by 

J°°(/) := a^x 1 A J , and J k (f ) := a !J i ! A J for a tJ G K. 

(i,j) i+j<k 
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The set of terms of / are defined by 

Terms(/) := {a^x 1 X J \ a^ 0} , 

i.e., all terms appearing in When Terms(/) 7 ^ 0, the leading term of / is 

LT(f) = max Terms(/), 

i.e., LT(f) G Terms(/) and p -< LT(f) for any p G Terms(/) \ {LT(f)}. The germ g is flat iff 
Terms(/) = 0. We define LT(f ) := 0, when / is a flat germ, i.e., J°°(/) = 0. The coefficient 
and monomial of the leading term are respectively called the leading coefficient ( LC ) and leading 
monomial ( LM ). For the case of / G K[x, A] and -<i ex , we may similarly define LT(f),LC(f) and 
LM(f). 

Now we present some known definitions, terminologies and theorems from (6j[lHIT5|l32] . These 
are suitably modified and generalized to fit in our purpose. 

Definition 2.3 (Remainder). A remainder of a germ / with respect to the set of germs G = 
{gx, ..., g m } C S and the local order -< is defined as a germ Rem(/, G, -<) G S so that 

( 1 ) Rem(/, G, -<) = f - Yn=\ Q.%9h for some ?i, ■ • •, Qm e S so that LT^gf) < LT(f). 

(2) No term of Rem(/, G, -<) is divisible by any of LT(gf) for i = 1,..., m. 

By replacing S with K[[x, A]] in Definition 12.31 the remainders in the ring of K[[x, A]] is readily 
defined; also see [ 6 | Pages 251-252], The same is true for the case of polynomial germ ring K[x, A] 
provided that the local order -< would be instead a monomial ordering like -<i ex - The remainders 
in the ring of SS is defined in m Page 170] using Mora normal form algorithm. Usually, the 
terminology of normal form is used rather than remainder. However, we choose remainder as the 
other may cause confusion with the normal forms of germs in singularity theory. The division here 
is related to Malgrange preparation theorem and Mather division theorem; see (41] Corollaries A. 6 .2 
and A.7. 2 , Theorem A.7.1]. 

When / is flat, its remainder with respect to any set of germs and local order is flat. The 
remainder is not necessarily unique even modulo flat germs; also see [ 6 j Pages 251-252]. In fact, the 
remainder is unique modulo flat germs when G is a standard basis (or Grobner basis for the case of 
K[x, A]), defined as follows. 

For an ideal / in 77, we define the leading term ideal LT(I)n by 

LT(I) n := (. LT(f) : / G /)*. 

Definition 2.4 (Standard basis). Let I be an ideal in 77 with a finite generating set {g ±,..., g m } C 
I. The set {rq,..., g m } is called a standard basis of / when 


LT(I) n =(LT( gi ),...,LT(g m )) n . 


(2.1) 
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The set {g \,..., g m } TZ is called a standard basis in TZ when it is a standard basis for the ideal 

{dli • • • i drain- 

Remark 2.5. (a) The set {gi,..., g m } C K[x, A] is called Grobner basis with respect to -<i ex 

when the condition ( 12.1 ft holds. 

(b) Any finite set of germs {gi | 1 < i < m} in S’ is a standard basis in S iff the set of formal 
power series {J°°(gi)} is a standard basis in K[[x, A]]. 

(c) Let 7 Z := K[x, A], K[[x, A]] or M and {g ±,..., g m } C / be a finite set. Then, 

LT(I)k — {LT{g l ),..., LT(g m ))iz 

implies that {gi ,..., g m } C / is a generating set for I. This has a simple argument as follows; 
also see pH Page 206]. For any / G I, 

r := Rem (/, {&| i = 1,..., m}, -< ) el. 

Thus, LT(r) G LT(I)n and LT(r ) must be factored by LT(gi) for some i. The latter implies 
that r = 0 , otherwise this contradicts with r being a remainder. 


Computation of standard basis uses the notion of S'-germs. Let f,ge7Z and ^ be a local 
order. Then, S-germ of / and g is defined by 

f LCM{LM(f),LM(g)) LCM ( L M(f),LM{g)) + P 

S(f,g)=< LT (f) * lt(9) 9 1 J U> J ^ g ’ T u > 

{ 0 for J°°(/) J°°(g) = 0. 

Here LCM stands for the least common multiple for a pair of monomials. 

Theorem 2.6. (Also see Hironaka Theorem on /S Page 252]) Let TZ = S or K[[x, A]], G = 
{gi ,..., g m } C TZ, 0 f eTZ and -< be a local order. Then, 

(a) Always f has a remainder with respect to -< and G. 

(b) The set G is a standard basis iff the remainder of f with respect to G and -< is unique modulo 
flat germs. 

(c) (Buchberger’s Criterion) The set G = {gi,...,g m } is a standard basis iff for all i,j, the 
expression Rem (S(gi, gj),G, -<) is a flat germ. 


(d) The set G is a standard basis 'iff Rem(/, G, -<) is flat for all f G {G)k. 
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Proof. When 1Z := K[[x, A]], the claim is Hironaka Theorem given on p, Page 252]; also see [35] . 
Therefore, we assume that TZ — $. By Hironaka Theorem, the remainder always exists in the ring 
of K[[x, A]], i.e., there exists a 

Rem(j“(/),{J“(9,)}£„^ ) £ R[M]] 

and Qi G K[[x, A]] so that 

m 

Rem(J~(/), -< ) = ~ 

i— 1 

with LT(J°°(/)) -< LT(gjJ°°((?*)). Further for 

P e Terms^Rem(j°°(/),{J 0O (^)}^i,^ )), 

p is not divisible by any of LT (J°° (gfl) for i = l,... pm. Our proposed division algorithm for germs 
in <§ follows 

m 

Rem(/, { gi}T=n-<) '■= f - 

i= 1 

= Rem(J“(/),{J“(ffi)}ti,-<) ( 2 . 2 ) 

m m 

+/ - j~(/) + T - 5] *9., 

2=1 2=1 

where g* G J°°(qi ) := qi- Note that the last line in (12.21) is merely a flat germ. Equation (12.21) 
implies that the division in S’ merely differs by flat functions from its analog in K[[x, A]]. The proof 
of part (a) is complete by Equation (12. 2 j) . 

Let r = f — Y1T=i QiSi satisfying (1) and (2) in Definition 12.31 Then, 

m m 

= j~(/) - J”toi9 ( ) = /“(/) - 

2=1 2=1 

= Rem 

Therefore, Rem(J°°(/), {J°°(g i )}”L 1 ,-<) is unique if and only if {J°°(gi)} is a standard basis in 
K[[x, A]]. This completes the proof by Remark 12.51 The proof of parts (c) and (d) are given by 
J°°{S{g il gj)) = S(J°°gi, J°°gj), Equation (12.2[) and Remark [2.51 a). □ 

In order to compute the standard basis for an ideal in K[[x, A]] (and also in S’), one needs 
to sequentially enlarge and update its set of generators by only adding the non-flat remainders 
(with respect to the updated generators and a given local order) of the S-germs of the generator 
pairs. This process is usually known as Buchberger algorithm and it terminates when the ascending 
ideals generated by the leading terms of the updated generators stops any further enlargement. The 
Buchberger algorithm is finitely terminated due to the fact that K[[x, A]] is a Noetherian ring. 
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Example 2.7. (a) The division of a polynomial by a set of polynomials in our division algo¬ 

rithm may involve formal power series; for example we have Rem(l, 1 — x, -< a iex) = 0 = 

1 “(E^o *")(!-*)■ 

(b) The remainder of a polynomial divided by a polynomial may give rise to an infinite formal 
power series. For instance let G := {xA — x 2 A 2 — x 4 }. Then, 

h := Rem(x 2 A, G, -< a iex) = x 5 + x 9 + 2x 13 + 5x 17 + 14a; 21 + ... 

= x 2 X — (x + x 2 X + x 5 + a: 3 A 2 + 2a; 6 A + a; 4 A 3 + ...) (xA — x 2 X 2 — x 4 ). 

Since the generator of an ideal with only one generator is always a standard basis, the re¬ 
mainder here is unique. 

(c) This example is to show that a finite set of polynomial ideal basis in S may lead to a standard 
basis that includes non-polynomial germs. Let G = {/ := xX — x 2 A 2 — x 4 , g := A — xA — xA 2 — 
x 3 }. Then, S(f,g) = x 2 A. It is easy to verify that S = {/, g, h} is a standard basis, where h 
is defined in part (b). 

The first and second Buchberger criteria are applied for efficient computation of standard basis 
in Singularity. However, we will not discuss them in this paper. 

Definition 2.8 (Reduced standard germ basis). Let G = {gi,...,g n } be a standard basis and 
LC(gi) — 1 for i — 1,..., n. When 

• LT(gi) \p for all p G Terms [g-j ), except for when p = LT(g t ) and i = j, 
the set G is called a reduced standard basis. 

Given a local order -<, [6j Theorem 2.1, Page 255] states that any ideal in K[[x, A]] has a unique 
reduced standard basis. 

Theorem 2.9. With respect to any local order -<, every finitely generated ideal ICS has a reduced 
standard basis. The standard basis is unique modulo flat germs. 

Proof. Put / = J°°(J) := ( J°°f | / G I). By [6, Theorem 2.1, Page 255], there exists a unique 
reduced standard basis {gi,... ,g m } for / where g % = J°°(Gi) with G* G I. Since LT(Gi ) = LT^gfl 
for i = 1 ,..., m, 

LT(I) K [[ Xi x]] = ( LT(gi ) : % — 1,..., m) K [[ X}X ]] = (LT(Gf) : i — 1,..., m) K [[x,x]]- 

{- LT(f) | f el} = { LT{f) | / G /} C (LT(Gi) : i = 1,..., m) K[[x , x]] 

C (LT(Gi) — m)g. 


Now we have 
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Hence, LT(I)g = (LT(Gi) : i = l,...,m) < g i and (LT(Gi) : i = 1 C I. Assuming that 

/ := ({fj}" =1 )g-, the set {G i } r [L l U {Rem(/j, {G i ]'^ =ll is a reduced standard basis, since 

Rem(/j, {GiYff 1; -<) is flat for any j due to part (c) in Theorem 12.61 The uniqueness follows [6] 
Theorem 2.1] and Remark 12.51 fab □ 

The following Theorem along with Equation (I2.2[) and Buchberger algorithm provide compu¬ 
tational guidelines on how to compute a reduced standard basis in S. 

Theorem 2.10. Let I = </i,..., f n )g. 

• There always exist germs (?i,..., g m G S for m < n so that 

I = (9u---,9rn)* (2-3) 

and LT(g.i) \p for all p G Terms^-), except for when i = j and p = LT^gf). 

• Furthermore, assume that {/,;}’L 1 is a standard basis in S. Then, the set (form < n) 

can be chosen so that it is a reduced standard basis. 

Proof. The proof here is similar to the idea on [6] Page 255] and gives us an algorithm for reduced 
standard basis computation. Without loss of generality assume that LT{ff) \ LT(f 3 ) for i ^ j\ 
otherwise, fj = Yhi^jTfi + r j f° r some r 3 , (p G S satisfying conditions in Definition 12.31 Thus, fj 
may be replaced by r 3 and this may cause m < n for when r 3 = 0. By Theorem 12.61 there exist 
qij G S for j — 1,..., n so that 

r i = Rem (/, - LT(fi), {f 3 }] =11 F)=f i - LT(ff) - ^ q, 3 f 3 . 

Now for each p G Terms(rj) we have LT(fj) \p. Define gi = LT(fi) + r*. This shows gt G / and the 
first claim follows. The second claim is valid due to LT(/j) = LT^gf). □ 

Example 2.11. Let / be a non-zero flat germ, and / = ( G)g where 

G := {g i := A — A exp (x), g 2 := x — s'm(x), g 3 := Xx + A 3 + A 2 ln(l + x) + /}. 

Since LT(G ) = {Ax, x 3 } and LT(I ) = (Ax, x 3 , A 3 )^, G is not a standard basis with respect to -< a iex- 
Dne to 

Rem(S(g 1 ,g 2 ),G, -C a i ex ) = 0, Rem(S(g 2 , g 3 ),G, -C a i ex ) = -x 2 / 

and Rem(S'((yfi, ^ 3 ), G, -< a i ex ) = A 3 + /, we add A 3 + / to the basis, i.e., SA := {^i, ^ 2 , ^ 3 , A 3 + /}. 
Now Si is a standard basis, yet it is not a reduced standard basis. Given Theorem 12.101 and its 
proof, 

LT{gf) | LT(g 3 ), Rem(. 93 , Si \ { 93 }, -< aiex ) = 0, S 2 := { 91 , g 2 , A 3 + /} and 

Rem (Ax + g x , S 2 , -< a i ex ) = Rem [g 2 - S 2 , -< alex ) = 0, 

we define 

S 3 := {Ax, x 3 , A 3 + /}. 

Now S 3 is a reduced standard basis. 
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2.2 Finite codimension 

Finite codimensional ideals demonstrate an important role in ideal presentation in Section [3] Let 
M-g (x, X)$, = (x, X)&, Mk[[x,\}\ — (x, A)a-[[ x ,a]]) and M.k[ X) \\ = (x, A)#^. 

Definition 2.12. An ideal / in S (or in PS, K[[x, A]]) is said to have a finite codimension when 
M.g k Q I (A4^ k Q I, M-K[[x,\] k C I) for k G N. Equivalently, I has a complement (vector) subspace 
in S with finite dimension. 

Now we compare the ideals in K[x, A] with those in K[[x, A]] and S. 

Theorem 2.13. Let I be a finite codimensional ideal in S. 

(1) The ideal I has a unique reduced standard polynomial germ basis. 

(2) Assume thatG := {c/j}” =1 is a standard basis for I and f 6 S. Then, the remainder is always a 
polynomial germ, i.e., Rem (f,G, -< a iex) € K[x,\\. In particular, we have deg(Rem (f,G, -< a iex 
)) < k when M/ +1 C /. 

Proof. By Theorem 12.91 we assume that G := {/)}( l =1 is a reduced standard basis for /. Assume that 
M.g k+1 C / for some k. Then, for any 0 < l < k + 1 there exists fieGso that LT(ff)\x k+1 ~ l X l . On 
the other hand for any p € Terms (fj — LT(fj)) and j < n, LT(ff) \p and hence, x k+l ~ l X l \p. The 
latter implies that J°°(ff) = J k+1 (fj ) for all j. Finally, we recall that any finite codimension ideal 
includes the set of all flat germs. This allows us to replace fj with gj := J k+1 (fj ) and S := {gj}f =l 
is a reduced standard basis. The second claim similarly follows from A4g k+1 C LT(I)g = ( LT(gf))g 
and the fact that no term of Rem(/, G, -< a iex) is divisible by LT(gf). □ 

Example 2.14. Let S = {./)}" =1 be a reduced standard basis whose generated ideal is of finite 
codimension. Then, Rem(/, S, -<! a i ex ) = / — J^difi £ K[x,X]. However, q j may not always be a 
polynomial germ. Consider the ordering A -< a i ex % and define 

I := (gi := 2A 3 — 3X 2 x + x 5 , g 2 : = — 3xX 2 + 5x 5 , g 3 := —3A 3 + 5x 4 A) . 

The ideal / is finite codimensional since A4g 6 C I. Here, G = {gi,g 2 , < 73 , §:r 5 — ^x 4 A} is a standard 
basis and 

S = |/i := xX 2 - y^x 4 A, f 2 := ~g 3 , f 3 := x 5 - j|x 4 A J 

is a reduced standard basis. Now we have 

x 3 A 3 = x 3 f 2 + + ^A x 2 qf 3 , 

where q : = ^°^ 0 (y§§)* 2 ; 2 \ Rem(x 3 A 3 , S, -< a i ex ) = 0 while q\ := j^x 5 q and q 2 '■= |A x 2 q are not 
polynomial germs. 
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The following theorem enables us to compute the standard basis for ideals in S’ through com¬ 
putations in the fractional germs. 

Theorem 2.15. Consider I = (p±,... ,p n )g, -< a iex and M.^ k C (p 1; ... ,p n )& C Further, suppose 
that {q \,..., q m } is a standard basis for (pi,... ,p n )@- Then, {qi,..., q m } is a standard basis for I. 

Proof. See Appendix 17.11 □ 

The following theorem is one of our main contributions in this paper and provides important 
alternatives and criteria for our computations in different circumstances including ideals with infinite 
codimension. 

Theorem 2.16. (1) Suppose that G {gi,...,g n } C Then, the ideal (gi,■ ■ ■ , g n )s has a 

finite codimension iff {gi, ■■■, g n )sg I s a finite codimension ideal. Assuming that (gi,, g n )& 
has a finite codimension, for nonnegative integers i and j, JAg l (Xffg C {gi, ■ ■ ■, g n )s iff 
M@(\ j )& C (g 1} ..., g n )&. 

(zl) Let AAk[x,x] — (gii • • • •> 9n) k[x, a] — A [x, A]. Then, (XI s ) k[x,x\ — (9h • • • i gn) k[x, a] iff 

C (g 1 ,..., g n )g. 

(3) For a finite codimension ideal (gi ,..., g n )$g in and an ideal I in S , 

I = (gi, ■ ■ ■ ,9n)g iff I G & = (gi,...,g n )<%. 

In particular, let /, J be two ideals in S’, I have a finite codimension and I fl FS = J fl 
Then, I = J. 

(4) For an ideal I in S’, the following three conditions are equivalent. 

(i) M / C I, (ii) M,% k cm&, (ill) M k[x , A] fc C / n K[x, A], 

(5) Let gi,..., g n G S’ and k ^ N, for k, N E N. Then, either of the conditions 

Myff C (J N g 1 ,..., J N g n )>% and M K [[x,\]} k c (J N gi ,..., J N g n ) k[[x,x\\, 
is equivalent to M.g k C (g 1 ,..., g n )g. 

(6) Consider the finite sequence of nonnegative integers fffii E N so that the sequence ff + l t is 
decreasing and h is increasing. Let (ff, f 2 ,..., f n )g be an ideal that is not necessarily of finite 
codimension and 

I := ^AV"(A% C ..., J N f n ) M . (2.4) 

Then, either of the following conditions 

— for each j < n, Terms (fj — J Pj ff) C A iggl for some pj < N. 
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— for each j < n, fj G and fj — J Pj fj G A iggl for some pj < N. 
implies that 

(2.5) 

Proof. Part (1). The if part is trivial. Thus, we assume that Aig k C (g 1: ... ,g n )g and prove that 
M.@ k C (gi ,..., g n )gg. Using Nakayama lemma [29, Lemma 5.3, Page 71], it is enough to verify 
that M@ k C (gi,..., g n )gg + M,% k+1 . Since Mg k C (gi,..., g n )g, for some a s G S we have 

n n n n n 

X i x k-i = J2a s g s = J k J2^gs = ^ jk ( J ^ a M=J k J2 b ^s = J2 bs 9s + P 

5=1 5=1 S= 1 5=1 5=1 

where a s G S, b s : = J k (a s ) and r G M.@ k+1 . For the second claim, we do not need Nakayama 
lemma. We assume that M.@ k C (g 1 ,..., g n )&. Similar to above, for i + j < k we have x l \ l+ ^~ l = 
J k ^ 1 (x l = J2s =i bs9s + r i where r G and h s G 

Part (2). Proof is similar to the proof of the second claim in part (1). 

Part (3). Assume that / = (< 71 ,... ,g n )g- Trivially, 

{dli • • • 1 9n)^ L! I n (< 7 l, • • • , 9n)&- 

Let / G / fl &. Thus, / = Yh a i9i f° r some a* G S. So, / = J k Y J k ( a i)9i + h for an h G M-K[x,\ k+1 ■ 
Since the left hand side belongs to (< 71 ,... ,g n )&, the proof of the if part is complete. Now assume 
that I = (gi ,..., g n )& and / G I. By part (1) and (gi ,..., g n )g C /, J k (f) G / for some k gN. 
Since J k (f ) G / fl 

jk U) = 22 aiSi ’ for ai e ^ 

On the other hand / — J k (f) G (< 71 ,..., g n )g. This completes the proof of part (3). Part (4) is 
trivial. 

Part (5). Now the assumption A4gg k C ( J N g \,..., J N g n )s implies that for some a* G and 
any % — 1 ,... k, 

n n n 

x k ~ l X l = J k 22 0 ,iJ N {gi) = J k 22 ai 9i = 22 aiSi + r ’ for k ^ N} 

i =1 i =1 2=1 

where r G A4,^ k+1 . This and Nakayama lemma conclude that M.g k C (< 7 ^... ,g n )g. The converse 
and rest of the proof use similar arguments. Note that Nakayama lemma given in [291 Lemma 5.3] 
is also true when S is replaced with ^ and K[[x, A]]. 

Part ( 6 ) is concluded by Nakayama lemma and the fact that x ki ~ m \ li+m = Y aiJ N fi = Y ctifi + 
r, where r G (Terms (fj - J N fj))<% C (Terms(/j - J Pj fj))<% C M#I. □ 

The hypothesis of part (1) in Theorem 12. 161 allows us to use fractional germs for computations 


while the condition in part (2) permits the use of Grobner basis. 
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Example 2.17. Part (1) from the previous theorem is not valid when £% is replaced with K[x, A]. 
For instance, consider the example given in [38] Table 1, III. 1 for k — 5] and [221 Page 77], 

I := (gi := x 5 + x 3 A + A 2 , g 2 : = 5x 5 + 3x 3 X, g 3 := 5x 4 X + 3x 2 X 2 )g. 

The ideal / has a finite codimension since A4# 6 C /; see [29, Page 77]. However, the ideal Iq : = 
(gi, g 2 , g 3 )k[x,x\ h as an infinite codimension in K[x, A]. The reason for this is as follows. The reduced 
Grobner basis of / with respect to A -<i ex x is given by 

G := {3125A 3 + 108A 4 ,18A 3 + 125A 2 x, 2a; 3 A + 5A 2 , 2x 5 - 3A 2 }. 

Further for any natural number n > 3, the remainder Rem(x n , G, -<i ex ) = cA 3 where c G K. 
Therefore, x n does not belong to Jo and Jo is an infinite codimensional ideal. 

3 Finite codimension ideal representation using intrinsic 
ideals 

An elegant approach for answering the ideal membership problem is to provide a good presentation 
for ideals. In this section we define intrinsic ideals and use it for such representation. This plays 
a central role in developing Singularity and is a prerequisite for most computations of objects 
presented in Section [4] 

Definition 3.1. When only the trivial reparametrization A(A) = A is allowed in Equation ( 11. 2 j) . 
the associated relation is called strongly equivalent relation ; see [29, Page 51]. Then, an ideal / 
is called intrinsic [29, Page 81] when / includes all strongly equivalent classes of its members. 

By [29l Proposition 7.1], every finite codimension ideal / in $ is intrinsic if and only if there 
exist nonnegative integers s, ?»,, Ui for i — 0, ..., s so that 

S 

I = A ni }^, (3.1) 

i=0 

71q = 0, and the sequence Hi is strictly increasing while mi + Ui is strictly decreasing. The conditions 
on mi and n t make the representation (13.11) unique. Equation (13.11) for intrinsic ideals gives a 
convenient answer for the ideal membership problem. The monomials x mi X 71i for % = 0,..., s are 
called intrinsic generators of the intrinsic ideal /. For non-intrinsic ideals or more generally for a 
vector space /, we define their intrinsic part; i.e., we denote Itr(J) for the largest intrinsic ideal 
contained in /. 

Lemma 3.2. For an intrinsic ideal I, there always exists a reduced monomial standard basis for 
I that includes its intrinsic generators. Let I be a finite codimension vector subspace of S’ or be 
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a finitely generated ideal. Then, there exist nonnegative integers m and n, a reduced monomial 
standard basis {fj}j =1 for Itr(/) and {.9*} mi C S so that 

I = Itr (/) + {{ 9i }Zi)n, (3.2) 

where 7 Z := S or TZ := Ii. Here, none of the terms in Terms(r/j) is divisible by f) , i.e., Itr(/) D 
(U™ 0 Terms (^)) = 0. 

Proof. Let the intrinsic ideal / be given by (13.ip and G be the set of its intrinsic generators. Then 
for i = 0,..., s, we consecutively update G by adding the monomials of type x a \ 13 to G where 
a + = mi + 7ii, > Hi, and x a X 13 is not divisible by the elements of G. This gives rise to a reduced 
standard monomial basis for I. For the second part, we divide the generators of / by the reduced 
standard monomial basis of Itr(J) and define the nonzero remainders as g*. □ 

A refinement of the decomposition (13.2p for finite codimension ideals is given as follows. We 
denote I 1 for the set of all monomials not in /, while stands for the vector space generated 

by / x . Now by [[29] Corollary 7.4], we have 

I = Itr(/) © (7 n (Itr(/) ± )• (3.3) 


3.1 Computation 

In this section we describe our suggestion on how to compute the intrinsic ideal representation (13.ip 
and next, the representation (13.31) for finite codimensional ideals. The first step here is to find a 
lower and upper bound for the maximum natural number k so that A4^ fc C / < S. 

For any ideal J in a ring 1Z we define 

Tuj ■ y —> y, Tu,j(f + J) '■= uf + J, (3.4) 

for u G {x, A}. Obviously for 

J ■= I and 1Z \= S, 

S // is a finite dimensional vector space and is a linear nilpotent map of nilpotent degree 
N u < k. Therefore, we have 

max{W, N\} < k < N x + N x . 

Now we explain how to derive N x and N x . Computing a (reduced) standard basis for /, say S = 
{gi,i = 1,.. .n}, we now introduce a vector space basis generator 

B = {w | w is a monomial, w ^ ( LT(gi ) | gt G *5)^}. (3.5) 

for (f/J; also see [T4] Pages 177-179 and Proposition 4] and [3J, Pages 128-129]. 
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Lemma 3.3. Let I be a finite codimension ideal in S and S be a standard basis for I. Then, B +1 
is a K-vector space basis for S/I. 

Proof. Since / is a finite codimensional ideal, is so LT(I)g. Thus, B is a finite set. Assume 
that there exists aiWi G / for some a* G K, where Wi G B. Since S' is a standard basis, 
LT(f)\LT(Yf™ =1 aiwf) for some / G S. However, LT(f^fiLi a i w i) £ { a i w i}iLi- This contradicts with 
(13.511 . unless a; = 0 for all i. This concludes that B is K -linearly independent. 

For any g e S we have g + I — Rem(g, S, -< a iex) + T 

Terms(Rem(g, S, -< a i ex )) n LT(I) S = 0, 

and Rem(r/, S', -< a i ex ) is a polynomial germ. This concludes the proof. □ 

Lemma [3.31 readily provides a matrix representation (and thus, its nilpotent degree) for (p u j. 
The natural number N, satisfying niax{iV x , N\} < N < N x + N\ — 1, is consecutively increased 
to obtain k = mo. In fact, the remainders of x l X N ~ l divided by the standard basis S conclude the 
result, thanks to Theorem 12.61 (part d). 

Next we look for the maximum values of m; and rij such that A4g mi (X ni )g C I. Our suggestion 
uses the concept of colon ideals. Let Ji and J 2 be two ideals in S. We define the colon ideal of J\ 
by J 2 (see m Definition 5, Page 194]) as 

J\ '■ J 2 = (/ £ $ | f d 2 0 J\)g. 

By an inductive procedure and repeating the above for the colon ideal / : (X n )g, we obtain m, 
and rii as desired. (Note that 1 < n, and mj + < m.j_ 1 + rij_ 1 < m 0 for any i.) Therefore, 

the only remaining challenge is to compute the maximum value m so that M.g m C I : ( X n )g 
or equivalently, X4g m (X n )g C I. The following lemma and its’ follow up comments facilitate the 
computations of a finite generating set for colon ideals and next, their (reduced) standard basis 
readily determines m. The procedure mentioned above (using Equation (13.4p and J I : ( X n )g ) 
provides suitable lower and upper bound for m. 

Lemma 3.4. Suppose that I is a finite codimensional ideal in S , J I D K[x,X] and J = 
(01, ■ • •, 9u)k[x,\}- Then, I = {g u ..., g n )g. Further, let J n (g)K[x,\] = {hi, • • •, h n ) K [ x ,x] for a mono¬ 
mial germ g G K[x, A]. Then, 

1 ■ {g)s \ . (3.6) 

\ 9 9 / s 

Proof. Let / G I. Since / has a Enite codimension, A4g k+1 C / for some k. Hence, J k (f ) G 
/ fl K[x, X] = J. Theorem 12.161 (part 3) completes the proof of the first part. 

For any u G / : {g)g, we have u(g)g C I. Thus, 


ug el n {g) g . 


(3.7) 
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For any pel D {g)g D K[x, A], p = fg = J degp deg9 (f)g for some f E S. Thereby, 

I n (g)g n K[x, A] = J n {g) K [x,\\ = {hi, ■■■, hn) K [ Xi x]. 

By the first part, we have 

7n {g)g = {hi,...,hn)g. 

This and Equation (IX7D imply that u G (^,..., *f)g and / : {g)g C (^,..., ^)g. 

Now let u G (y,..., y)^. Hence, aug G {hi ,..., h n )g for any a G S. On the other hand, 
{hi ,..., h n )g C / shows that aug G / and eventually u G / : (g)^. □ 

Given Lemma 13.41 we would merely need to compute J D (. 9 ) ft'fx.A] f° r the monomials g := A n . 
Note that the above lemma is a slightly modified version of [21 Theorem 11, Page 196] and [T5l 
Theorem 5.5, Page 185]. We simply use J 0 {g) k[x,\] rather than using the naive and alternative 
choice / D {g)g- The reason for this is that the computation of J 0 {g) k[x,\] using the lexicographic 
ordering is efficient and classic (described below) while such an approach for I D {g)g does not seem 
easy to work with standard basis with local orderings. 

Finally, we recall the classical approach on the intersection computation of two ideals Ji and 
J 2 in the polynomial germ ring K[x, A]; see [2, Page 187] for more details. Let Ji := (/) | i = 
ni) K [x,\] and J 2 := {g t \i = 1 ,..., n 2 ) K [x,\}- We define 

J 3 ■= (tfi, (1 - t)gj | for i = 1,..., m and j = 1,... n 2 ) K[x X t] . (3.8) 

Thus, J 3 D K[x, A] = J\ fl J 2 . We may compute a Grobner basis for J 3 with respect to -q ex (where 
X , A “^lex t") 

G 3 := {pi(x, A), qj{x, A, t) | for i — 1,..., n and j = 1,... m}, 

here Pi represents those basis elements independent of t while qj represents those explicitly depending 
on t. Since {LT{G 3 ))K[x,\,t\ = LT(J 3 ), we may conclude that 

G := {pi{x, A) | for i — 1,..., n} (3.9) 

is a Grobner basis for J\ fl J 2 . The reason is as follows. For any u G LT(Ji fl J 2 ) C LT(J 3 ), either 
LT{pi ) or LT{qj ) for some i or j must divide u. Since t divides the leading term of qj for any j, 
u G {LT{G)) K [x,x]. 

Given Equation (13.31) . we merely need to use the generators of /, a standard basis for Itr(J) 
and Theorem 12.61 (part d) in order to obtain / fl (Itr(/) _L )x- This completes the required approach 
for ideal representation of any finite codimension ideal in S. 

Remark 3.5. Applications of the above procedure is when / is given by a finite set of generators 
{di}Z= l- Thus, it is essentially useful (when it is possible) to instead use a truncated Taylor expansion. 
This is justified when the hypothesis of Theorem 12.161 (part 5) is satisfied. Parts (1, 3) in Theorem 
12.161 indicate that for finite codimension ideals in S, we may simply replace 7 Z := S with either 
3% or K[[x, A]]. Part (6) in Theorem 12. 161 provides an important computational criterion for ideals 
with infinite codimension in S. 
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4 Computations of objects in singularity theory 

In this section we recall the algebraic tools and present our suggested approaches that are needed 
for computation of normal forms, universal unfolding, and persistent bifurcation diagram classi¬ 
fication. We acknowledge frequent fruitful discussions with Professor A. Hashemi. The second 
author acknowledges B. M.-Alizadeh’s helps, ideas and lessons in the early stages of this project. 
They were essentially helpful in a fast learning of the concepts from algebraic geometry and their 
programming. 

4.1 Normal form 

Given a singular germ g e S’, from m Pages 88-89] we recall the intrinsic ideals &(g) and S^{g) 
by 

&(g) ■= Itr {(xg, a g, X 2 g x , Xg x ) e ) , 

and 

S'(g):=Z (mi , ni )M, m *( X n %, (4.1) 

where ( 0 , 0 ) ^ 0 , and there would not exist nonnegative integers p and q such that 

-§^-§^{g) (0, 0) 7 ^ 0, q < Hi, p + q < m* + n*. The extra restrictions on m, and rii here make the 
presentation (14. ip unique as of those in Equation (13.ip . 

m Proposition 8 . 6 ] indicates that for any germ p G AA, g±P is contact-equivalent to g. Terms 
in &(g) are called high order terms ; see [291 Page 89]. Therefore, we may stay contact-equivalent 
to g by removing all terms in Terms(< 7 ) D ^(g) from Taylor expansion of g. Therefore, &(g) has a 
finite codimension if and only if g is finitely determined, i.e., g is contact-equivalent to a polynomial 
germ. 

[ 29] Theorems 8.3 and 8.4] state that for any term x rn X n in ^(g) 2 -, we must have J-yrJyrr (g) (0, 0) 
= 0, while -§^-§^{ 9 ) (0,0) 7 ^ 0 for any intrinsic generator x m X n in S^{g). Now we recall the inter¬ 
mediate order terms as terms belonging to 

A := f^ > (g) ± \ (fA ± (g) {x mi A ni | x mi X ni is an intrinsic generator for y(g)}^J. 

4.1.1 Normal form computation 

Now we are ready to provide an algorithm for computing normal form of a given germ g. Using the 
procedure given in Subsection 13.11 we may compute &(g) and remove all terms in Terms( 77 ) f| ^(g) 
from Taylor expansion of g to obtain a more simplified contact-equivalent germ, say /. Now it only 
remains to eliminate intermediate order terms A from / as many as possible. Then, this gives rise 
to its normal form. 

When A is empty, / will be called normal form of g. Otherwise, we may use suitable polynomial 
change of variable X(x,X) and positive polynomial germ S(x, A) to eliminate some intermediate 
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order terms from /. For example, we may replace x by ax + b\ in / where a > 0 and b are 
arbitrary constant coefficients. This gives rise to a system of linear equations and a maximal 
solvable subsystem leads to the further elimination of negligible terms in / and thus, the normal 
form computation of g. This approach needs to be systematically adopted along with standard basis 
computations for high order terms in multi-state variable cases. This is because complete algebraic 
characterization for high order terms is not yet known for many multi-state variable cases. This 
will be addressed in details in our upcoming result for multi-state variable cases. 

4.2 Universal unfolding 

In this section we recall the algebraic formulation needed for computation of the universal unfolding 
of a singular germ g. 

The tangent space T(g) is defined by 

T{g) ■■= (g, xg x , \g x )g © K{g x , g x , A g x ,..., A £ g x }, 

for sufficiently large l. For a finite codimension g, there exists a natural number t so that A £ g x ^ 
(g,xg x , A g x )f and A l g x G {g,xg x , A g x )g for l > t see [29] Page 127]. 

Next, a universal unfolding of g is defined by 

k 

G 0, A, /i) := g(x, A) + ^ a iPi( x i A), (4.2) 

2—1 

where a := (cq,... ,«*,), Pi-s form a basis for a complement space of T(g). Thus, we may choose 
pi G T(g) 1 -. The number k is called codimension of T(g) or equivalently codimension of g. 

The tangent space T(g) is computed via Lemma 13.21 arid the procedure given in Section l3d~1 
Since Equation (13.2p is a direct summation here, a complement space basis can be obtained through 
Equation (13.5jl and a linear algebra computation. 

4.3 Persistent bifurcation diagram classification 

Given our description in Section [T] persistent bifurcation diagram classification is complete by 
simplifying the defining equations of transition set £ and then, choosing one parameter from each 
connected component of £ c ; see [29], Page 140]. The latter is enabled in Singularity by using the 
Maple package RegularChains. Let g be a singular germ of finite codimension and G(x, A, a), a G 
K fc , be a universal unfolding of the germ g. We recall that the transition set £ := SB U U 
where 

:= {a G | G = G x = G x = 0 at (x, A, a) for some (x, A) Glx M}, 

:= {a G | G = G x = G xx = 0 at (x, A, a) for some (x, A) Glx M}, 

:= {a G | G = G x = 0 at (x i: A, a ) for i = 1,2 and X\ ^ a: 2 }. 


(4.3) 
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Possible reduction of variables in the defining equations (in particular removing x and A from the 
equations) in (j4.3j) is desirable. In order to achieve this goal, we let / C K[x,X,a] and J C ^ x ,\,a 
be the ideals generated by polynomial defining equations (14.3ft for either of SS and . Here <e x ,\,a 
stands for all germs of smooth functions of x, A, and a. For the case of 3>, a new variable ( is 
introduced and the quadratic germ 1 — <C(xi — x 2 ) is also added to the generators of the ideal / 
and J. This is due to the fact that x\ ^ x 2 . Next, we compute the Grobner basis G for I in 
K[x, A, a], a G M fc , with respect to aq -q ex A -q ex x for j = 1 Thus by [UJ Theorem 2], 

G D K[a\ is a Grobner basis for / D K[a\. Hence, G D K[a] C Additional restrictions on the 

transition sets are also obtained by the other elements of Grobner basis G. 

5 Main features of Singularity 

The reader are referred to our user-guide [21] for a list of all functions, their capabilities, options 
and comprehensive information on how to work with Singularity. In this section we merely 
describe the main features of Singularity. The main functions (not all) are given in Tables [H 
and [2] To illustrate how these functions work, let g(x, A) := x 5 + Ax + A 2 . Then, AlgObjects(g) 
returns & := M 6 + M 2 ( A) + (A 2 ), RT := M 5 + M (A), T := M 5 + M (A) + K{x + 2A,x 4 + ±A}, 
ff/T := K{ l,A,x 2 ,x 3 }, y := M 5 + M(X), := K{ 1, A, x, x 2 , x 3 , x 4 }, and intrinsic generators 

of 5? are x 5 and xA. For an example of computations in infinite codimensional ideals, consider 
restricted tangent space of g(x, A) = A 3 sin(x). RT(g, InfCodim) generates the restricted tangent 
space of g as A4(A 3 ). 


Table 1: Singularity’s main functions in singularity theory. 


Function 

Description 

Verify 

information on truncation degree and computational ring. 

AlgObjects 

2P, RT, T, $/T, y intrinsic generators of 5? . 

Normalform 

normal forms of a given germ. 

UniversalUnfolding 

universal unfoldings of a given germ. 

RecognitionProblem 

for normal forms and universal unfoldings. 

TransitionSet 

transition set are computed and plotted or animated. 

PersistentDiagram 

plots or animates all persistent bifurcation diagrams. 

Transformation 

estimates S and X relating two contact-equivalent germs. 

Intrinsic 

intrinsic part of a given ideal or vector subspace of §. 


Verify(g) checks a germ g for its bifurcation analysis while Verify(G, Ideal, Vars) checks 
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germ generators G of an ideal {G)g for divisions or standard basis computations. In either case, 
it suggests suitable computational rings and a truncation degree and it verifies that their use does 
not lead to error according to Theorem 12.161 

Morraalf orra(g) gives rise to the normal form x 5 + Ax while the universal unfolding 

G i ! — x^ T Ax T an T a? 2 A T 03 X T o 4 X^ (5.1) 


is derived by UniversalUnf olding(g) . The command TransitionSet (G_l) derives the transition 

set associated with G\ as 38 := {(an, a 2 , 03 , a 4 ) | an = a| + a^an — a|af 3 }, and 33 given in 
Appendix 17721 



Figure 2: Persistent bifurcation diagrams for parametric singularity G 2 in (15.21) . 


Now we consider the parametric germ 

G 2 = x T Ax T Q!i T a? 2 A ~\~ o 3 X (5.2) 

which is a universal unfolding for codimension three germ x 4 + Ax. TransitionSet (G_2) gives rise 
to 3$ := {(an, a 2 , ol 3 ) | o\ + afa 3 + an = 0}, rtf’ := {(an, a 2 , 0 ^ 3 ) 1 1280^3 + 3 a\ + 72ana| + 432a 2 = 
0, «3 < 0}, and 33 := {( 01 , 02 , 03 ) \ o\ — 4a 1 = 0}. A list of persistent bifurcation diagrams is 
generated by the command PersistentDiagram(G 2 , plot, ShortList, [x, A, an, a 2 , 03 ]); Figure [2] 
is some inequivalent diagrams chosen from this list. 

Let 


G := {g\ := sin(A 7 + x) + exp(x 4 ) — x — 1 — A 9 , g 2 := x 5 — A 2 , #3 := cos(x 6 ) — A — 1}, 
Vars := [x, A] and k be a sufficiently large truncation degree. Then, using the commands 
StandardBasis(G, Vars, k, Fractional), StandardBasis(G, Vars, k, Formal), 
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Table 2: Main functions of Singularity associated with the local rings K[[x, and <§. 


Function 

Description 

StandardBasis 

standard basis in either of the rings. 

Colonldeal 

computes the colon ideal given in Equation (13.61). 

Division 

remainder of a germ g divided by a set G. 

MultMatrix 

computes the matrix <p u ,j in Equations (13.41). 

Normalset 

finds a basis for complement space of an ideal /. 


and StandardBasis(G, Vars, k, SmoothGerms), we obtain the Standard basis of G in the rings of 
fractional maps, formal power series and smooth germs S’. In this example, the standard basis in 
either of these rings is given by {A,x 3 }. Note that the options ForraalSeries and Fractional in 
StandardBasis use a truncated Taylor series expansion of germs (for non-polynomial and non- 
fractional germs) and are adopted according to Theorem 12.161 

6 Future works 

Our future and in progress projects are about multi-state dimensional singularities, singularities 
depending on symbolic coefficients and the study of moduli, equivariant singularities with compact 
symmetries, and an integration of Singularity with the normal form analysis of singular ordinary 
differential equations. As our research progresses, we will update Singularity and its user guide 

m- 

In this paper, we have considered scalar singularities. However, multi-state dimensional singu¬ 
larities frequently occur in applications and an enhanced capability of Singularity for their study 
is an important contribution. Such enhancement enforces implementation of extended concepts and 
tools for module structures; see um\- 

The study of local zero structures for parametric smooth maps with symbolic coefficients in the 
vicinity of their critical points is important for applications. Here, we refer to parameters as small 
control parameters or perturbations while symbolic coefficients refer to those coefficients assumed to 
be constant. Further, the complex dynamics associated with singular germs with positive modality 
requires certain techniques from algebraic geometry so that they can work with symbolic coefficients 
and generate their associated classifications. This requires some extensions of standard basis; in our 
opinion, extensions similar to the notions of comprehensive Grobner basis and Grobner system in 
the context of local rings of germs; see [121133113411361141)] . This is in progress and will be addressed 
in details as a separate work. 

The possible enhancement of Singularity for equivariant singularities with compact synmie- 
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tries seems a challenging task but it is indeed a worthwhile contribution. We believe this should be 
closely related to a local germ ring version of SAGBI basis; e.g., see [37]. 

The contact equivalence relation is compatible with the equivalence relation arising from local 
invertible changes of variables and time rescaling in differential equations around an equilibrium; 
see [22b2?ll44l - 143] . Thus, Singularity can be integrated with the normal form computations of 
differential equations. The aim is to locate the local bifurcation of singular differential equations 
in terms of the coefficients and parameters of the original systems. This contributes to the local 
bifurcation control and to design efficient controllers for uncontrollable engineering problems; see 

W- 
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7 APPENDIX 

7.1 Proof of Theorem 12.15 

Let g G LTg(I) = ( LT(f ) | / G I)g. Since Q (pi, ■ ■ ■ ,p n )^, part (1) in Theorem 12.161 implies 
that M.$ k C / and M.g k C LTg(I). Hence, without loss of generality we may assume that g = J k {g). 
Thus, 

9 = Y, ^ LT Ui) = J k Y <D LT (fi) = Y J k (*i)LT{fi) + r, 

where a* G S’, fi G I and r G A4# fc+1 . Now we only need to prove that LT(fi) G ( LT(qj),j = 
1,... m)g when deg(LT(/))) < k. 

Since /* G I = ( qj,j = 1,... m) s , fi = Jf for some b ljt G S’. Thus, 

j\n) = J k E Jt{ - b >s)<n = E + s 

where s G N4^ k+1 . Therefore, J k (fi) G ( qj,j = 1 .rri)&. On the other hand, 


LT(fi) — LT(J k (fi)) G LT 3l {q 1 ,...,q m ) 3l ={LT{q 1 ),...,LT{q m )) a 

C (LT(qi ),..., LT(q m ))g. 


This completes the proof. 
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7.2 Parameter spaces J4? and Q) for G\ given in Equation (15.11) 

:= {(«i, a 2 , 013 , cnt) | 0 = 3375a2 a I a 4 + 10125aia2 tt 4 + 16875«2 a 3 — 675«2 a 3 a 4 + 288a 2 a|a| + 
67500aia| a 3 a 4 — 900aia 2 «3«4 + &Q^oiia\ + 1080a 2 «3«4 — 32a|«4 + 45000a; : j ! a 2 a;4 + 13500ai« 2 a| + 
3300 aia 3 a|— 108 af+ 30000 a 4 Q; 3 a 4 + 50000 a 4 } and S> := {(«i, a 2 , a 3 , a 4 ) | 2000a2a4+111000o;2a3a4— 
6400a2a 3 tt4 + 64a 2 «I + 28000aia2 a 4 — 16875a2 a 3 — 43200a2 Q; 3 tt 4 + 5472a 2 «3a4 — 128a 3 a4 + 
45000aia2a3Q;4—69600aia 2 a3a|+256aia|+34020a 2 Q;|Q;4—1728a3a4+130000a4a 2 a|—81000aia 2 af+ 
43200ct!ia|a| — 5832af — 180000a3a 3 a4 + 200000af = 0}. 



